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Charge-stabilized colloidal spheres dispersed in weak 1:1 electrolytes are supposed to repel each 
other. Consequently, experimental evidence for anomalous long-ranged like-charged attractions 
induced by geometric confinement inspired a burst of activity. This has largely subsided because of 
nagging doubts regarding the experiments' reliability and interpretation. We describe a new class of 
thermodynamically self-consistent colloidal interaction measurements that confirm the appearance 
of pairwise attractions among colloidal spheres confined by one or two bounding walls. In addition 
to supporting previous claims for this as-yet unexplained effect, these measurements also cast new 
light on its mechanism. 
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INTRODUCTION 

A long-lived controversy was ignited twenty years ago 
by the suggestion Q, Q that similarly charged colloidal 
spheres need not repel each other as predicted by Poisson- 
Boltzmann mean field theory 0, , but rather might ex- 
perience a long-ranged attraction for each other under 
some circumstances. Interest in this problem deepened 
when direct measurements of colloidal interactions re- 
vealed just such attractions in micrometer-scale colloid 
in aqueous dispersions at extremely low ionic strength 
Subsequent measurements demonstrated that such 
anomalous like-charge attractions are only evident among 
spheres confined by nearby charged surfaces, and not 
otherwise jf| 0, ■ This observation effectively refuted 
the originally proposed mechanism for like-charge col- 
loidal attractions 0,0, and other mean-field mechanisms 
were excluded soon thereafter on theoretical grounds 
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When the search for more sophisticated attraction- 
generating mechanisms subsequently failed to reach con- 
sensus, the experimental evidence came under renewed 
critical scrutiny. Measurements of long-ranged attrac- 
tions performed with optical tweezers near a single 
charged wall 0] were demonstrated to have been sen- 
sitive to a previously unsuspected kinematic coupling 
mechanism llRL 
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Suspicion thus was cast on all 
interaction measurements based o n op tical tweezer ma- 
nipulation in confined geometries |lS| |. Complementary 
interaction measurements performed on colloidal disper- 
sions in equilibrium are immune to kinematic artifacts 
H El IH El E2- However, they obtain pair po- 
tentials by inverting measured pair correlation functions, 
a process involving various poorly controlled approxima- 
tions. It is conceivable that these methods could mis- 
interpret oscillatory many-body correlations observed as 
attractive or even oscillatory pair interactions j^. In- 
deed, when particular care was taken to avoid such ar- 



tifacts in measurements in a carefully prepared model 
system, no sign of anomalous attractions was seen |l6l . 
These observations raise a disturbing question: could the 
entire case for confinement-induced like-charge attrac- 
tions be based on experimental artifacts? 

This article describes a new series of equilibrium col- 
loidal interaction measurements featuring novel tests for 
thermodynamic self-consistency. These measurements 
explicitly address all of the aforementioned sources of 
experimental error and yield equilibrium pair potentials 
(but only where appropriate!) with quantitative error 
estimates. Their results confirm that confinement by 
one or two nearby glass walls induces long-range equilib- 
rium attractions between nearby pairs of charged spheres. 
Confinement-induced attractions appear both among the 
highly charged polystyrene sulfate spheres that were the 
subject of the original round of anomalous observations, 
and also between more weakly charged silica. Trends 
observed with variations in confinement and electrolyte 
concentration shed new light on the attractions' origin, 
suggesting a role for nonmonotonic correlations in the 
distribution of simple ions near charged surfaces. 



THE STRUCTURE OF COLLOIDAL 
MONOLAYERS 



Our colloidal interaction measurements follow the gen- 
eral approach pioneered by Kepler and Fraden [j| and 
Vondermassen et al. |l9| . in which digital video mi- 
croscopy is used to measure the distribution of spheres 
in a dispersion at equilibrium. Figure ^ shows our imple- 
mentation schematically. An aqueous charge-stabilized 
dispersion fills a hermetically sealed slit pore between a 
glass microscope slide and a coverslip. The confined dis- 
persion is allowed to equilibrate with reservoirs of mixed- 
bed ion exchange resin to a base concentration of roughly 
1 fjM. Controlling the pressure of a buffer gas in these 
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FIG. 1: Measuring the structure of colloidal monolayers. 

reservoirs also permits the spacing H between the walls 
to be adjusted and maintained constant over the course 
of an hour |8j. Residual contaminant ions are believed 
to consist of sodium ion leached from the glass, and car- 
bonate infiltrating from the atmosphere, both of which 
are monovalent. The glass surfaces develop large nega- 
tive charge densities that repel negatively charged 
colloidal spheres and prevent them from sticking under 
the influence of van der Waals attraction. Depending on 
the resulting balance of forces on the spheres, the disper- 
sion can be confined to a monolayer at height h above 
the lower surface. 

Spheres larger than a few hundred nanometers in di- 
ameter are readily imaged by conventional bright-held 
microscopy. A detail from a typical video micrograph 
of (j = 1.58 /zm diameter silica spheres appears in fig- 
ureQ The spheres' centers can be tracked with standard 
techniques of digital video microscopy Q, with accura- 
cies approaching Ax = Ay = 30 nm being achieved for 
these particles [23l Eij . The plot in figure Q shows the 
trajectory of a single sphere over one minute from the 
region indicated by the box overlaid on the micrograph. 

The in-plane positions, r,-(t), of spheres labeled by j 
in snapshots obtained over time t can be compiled into 
the time-dependent particle density 



N(t) 

p(r,t) = ^(r 



(1) 



The rest of our results are extracted from p(r, t). 

For example, individual trajectories can be analyzed 
with the Einstein-Smoluchowsky relation 



P(Sk\t) = exp 



(5 k - Vktf 
2D k t 



(2) 



which describes the probability of finding particles dis- 
placed by distance Sk = (*j(t) — i"j(0)) fe along the fc-th 
coordinate after time t. Fitting to equation yields 
the particles' diffusion coefficients D k and mean drift 
velocities Vk- For an equilibrated isotropic system, we 
expect identical diffusion coefficients in orthogonal direc- 
tions and no overall drift. These conditions are met for 



all of the data sets presented below, with maximum drift 
speeds below 0.3 /im/sec and typical speeds far smaller. 

Provided care is taken to account for the finite field 
of view and the varying number N(t) of particles within 
it 0, |U [2^ , p(r, t) can be summarized with the radial 
distribution function 



9(r) 



1 /p (r-r',t) p(r',t) 
A(r) 



(3) 



where the angle brackets indicate an average over the 
field of view, over angles, and over time, and where n = 
N/A is the areal density of N = (N(t)) particles in area 
A, and A(r) is the area within the field of view over which 
pairs separated by r might be found. 



LIQUID STRUCTURE INVERSION 



The Boltzmann formula, 

g(r) = exp(-/3w(r)), 



(4) 



relates the radial distribution function for an isotropic 
system in equilibrium to the potential of mean force w(r) 
associated with its structure. Here, (i^ 1 = ksT is the 
thermal energy scale at absolute temperature T. The po- 
tential of mean force can be identified with the system's 
underlying pair potential only in the limit of infinite di- 
lution, 



u(r) 



lim w(r). 

n— >0 



(5) 



At higher densities, simple crowding can induce layer- 
ing, and thus oscillatory correlations, even in a system 
whose pair interactions are monotonically repulsive. In- 
terpreting the effective inter-colloid interaction is still 
more problematic. The spheres' dynamics reflect not 
only their direct Coulomb repulsions, but also the in- 
fluence of a sea of atomic scale simple ions, whose distri- 
bution also depends on the spheres' comparatively enor- 
mous charges and excluded volumes. The effective in- 
teraction between two spheres reflects a thermodynamic 
average over the simple ions' degrees of freedom. This 
almost certainly will depend on the distribution of other 
spheres at higher sphere concentrations. Under such cir- 
cumstances, the effective pair potential would not be well 
defined. At lower concentrations, however, the disper- 
sion's free energy can be described as a superposition of 
pairwise interactions. 

For all of these reasons, nonmonotonic dependence of 
(3w(r) = — ln<7(r) on separation r need not signal the 
onset of attractive interactions. Particularly in systems 
with long-ranged repulsive interactions, care must be 
taken to correct for many-body correlations. Unfortu- 
nately, no exact relationship is known between u(r) and 
w(r) at finite concentrations, even if the functional form 
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of u(r) is available. Instead, two strategies for invert- 
ing g{r) have emerged, one involving molecular dynam- 
ics or Monte Carlo simulations to refine trial pair poten- 
tials 0, [2(| , and another exploiting results from liquid 
structure theory to correct for many-body correlations 
[iH Eof . The results from either approach may be iden- 
tified with the underlying pair potential thanks to Hen- 
derson's uniqueness theorem |27j . 

We will avail ourselves of the Ornstein-Zernicke liquid 
structure formalism to invert a(r) l28j . building upon 
the pioneering work of reference |20j|. When applied 
to the spheres in a colloidal dispersion, the Ornstein- 
Zernicke equation describes how effective interactions 
among neighboring spheres give rise to structural cor- 
relations. In principle, it describes a hierarchy of N- 
body correlations emerging from pairwise interactions. 
Truncating the hierarchy yields analytically tractable ap- 
proximations, whose predictions are increasingly accu- 
rate at lower densities. Two of these approximations, 
the hypernetted chain (HNC) and Percus-Yevick (PY) 
equations have been found to accurately describe the 
structure emerging from computer simulations of systems 
with long (HNC) and short-range (PY) interactions. For 
two-dimensional systems, these are most conveniently ex- 
pressed as 

8u(r) = Bw(r) + { n/(r) (HNC) (6) 



where the convolution integral 



I(r) = / \g{r') - 1 - nl(r)] 



1] dV (7) 



can be solved iteratively, starting with J(r) = [29|. 
Evaluating I{r) directly rather than with numerical 
Fourier transforms minimizes the sensitivity of u{r) to 
noise in g(r). This implementation has been shown to 
be both accurate and effective in previous related studies 
EII2II2I. 



INTERACTIONS AND THE DLVO THEORY 

Figure shows typical results for pair potentials ob- 
tained from measured radial distribution functions with 
both the HNC and PY approximations. The data plot- 
ted as circles in figure Eta) were obtained for silica 
spheres a — 1.58 /im in diameter in slit pore of heights 
H = 195 /im and H = 9 /im. Silica's density is twice 
that of water, and these spheres sediment into a mono- 
layer with their centers at h — 0.9 /im above the lower 
glass wall, with out-of-plane excursions estimated 
to be no greater than Sh = 0.1 /im. This system was 
originally proposed as a model for studyi ng a ttractions 
mediated by a single wall in equilibrium |l6(. Indeed, 
the data obtained for a confined monolayer at H = 9 /im 



exhibit a strong and long-ranged attraction [23J. The 
pair potential measured at H = 195 /im, however, is 
monotonically repulsive 0, |2j| . This observation raises 
substantial questions regarding the nature of the more 
distant wall's influence. 




FIG. 2: Measured interactions in confined monolayers of (a) 
silica spheres a = 1.58 /im in diameter and (b) polystyrene 
spheres a — 0.652 /im in diameter. The silica spheres are 
sedimented into a monolayer at height h = 0.9 /im above the 
lower wall. The two data sets were obtained at areal density 
na 2 = 0.0654 for H = 9 /im and no 2 = 0.0797 for H = 
195 /im. The polystyrene spheres, by contrast, are confined 
to the midplane between parallel glass walls separated by H = 



1.3 /tm, with na 



0.056. 



The purely repulsive potential is described very well 
by the screened-Coulomb form predicted by the clas- 
sic Derjaguin-Landau-Verwey-Overbeek linearized mean- 
field model for colloidal electrostatic interactions 0, [1 : 



Pu{r) = Z 2 Xi 



exp(«;a)\ exp(— nr) 



1 



(8) 



Here, Z is the effective valence of a sphere of radius 
a = a/2, Xb = /3eo/(47re) is the Bjerrum length for a 
medium of dielectric constant e at temperature T, where 
eo is the elementary charge, and k" 1 is the Debye-Huckcl 
screening length given by k 2 = 47tAb71o in an electrolyte 
with a concentration n of monovalent ions. Fitting to 
the H = 200 /im data in figure|2Ia) yields a charge num- 
ber Z — 6500 ± 1000, in good agreement with predic- 
tions of charge renormalization theory > an( A screening 
length kT 1 = 180±10 nm consistent with the system's es- 
timated micromolar ionic strength. Comparable results 
are obtained for monolayers at areal densities ranging 
from na 2 = 0.04 to na 2 = 0.10, suggesting that the result 
is independent of density, and that the liquid structure 
inversion correctly accounts for many-body correlations 
in this concentration range. All other results reported 
here were obtained under comparable conditions. 

The observation of DLVO-like repulsions in a weakly 
confined silica monolayer is consistent with previous re- 
ports on this system |l6j. It also demonstrates that our 
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methods do not necessarily yield nonmonotonic poten- 
tials in this range of experimental conditions. When 
viewed in this light, the appearance of an attractive min- 
imum in the pair potential for the more tightly confined 
but otherwise identical monolayer at H = 9 /im seems 
more credible that it otherwise might [23|,|3- The obser- 
vation of attractions in silica colloid breaks the monopoly 
on anomalous attractions held by the substantially more 
highly charged polystyrene sulfate spheres used in previ- 
ous studies]! B, SEl 

Such indirect verification does not make the result any 
less surprising, however. The potential's minimum is 
roughly 0.3 fc^T deep at a center-to-center separation of 
r = 2.4 /im. The interaction's attractive component thus 
is substantially longer ranged that the core electrostatic 
repulsion and measurably influences colloidal dynamics 
a distances extending to several screening lengths. This 
greatly exceeds the range of like-charge macromolecular 
attractions ascribed to polyvalent counterions, countc- 
rion correlations, or fluctuations in the counterion distri- 
bution. Still more puzzling is that a wall separated from 
the monolayer by nearly 8 can qualitatively transform 
the spheres' apparent pair potential. 

Comparably strong and long-ranged attractions are ev- 
ident in the data plotted in figure |5Jb), which were ob- 
tained for polystyrene spheres a — 0.652 fjxa in diame- 
ter confined to the midplane between glass walls sepa- 
rated by H = 1.3 fj.rn. This is consistent with all pre- 
vious observations of like-charge attractions in confined 
polystyrene 0, H , including those involving optical 
tweezers ||. 

As an additional reliability check, results for the 
polystyrene data are plotted using both the HNC and 
PY approximations. Their quantitative agreement sug- 
gests that the monolayer's areal density is low enough 
for the liquid structure formalism to account accurately 
for many-body correlations in g(r). Indeed, there is little 
difference between w(r) and u(r) for this data set. We 
calculate the difference A«i(r) between the HNC and PY 
approximations for each data set and add it in quadra- 
ture to other sources of uncertainty to estimate errors in 
the reported u(r). 

By far the largest source of error results from experi- 
mental uncertainties in g(r). These, in turn, result from 
errors in measuring particle position and from counting 
statistics. Assessing the latter turns out to be somewhat 
subtle and establishes the lowest practical areal density 
n at which a reliable measurement can be made. 

The subtlety hinges on the following question: How 
many snapshots are required to ascertain whether or not 
the particles interact at all? In other words, how many 
pairs would we expect to see at the center-to-ccnter sepa- 
ration r in a non-interacting system? Given a spatial res- 
olution dr for binning particle separations into the radial 
distribution function, this number is 2irn 2 Ardr. Typi- 
cally, the number N = nA = 0(100) of particles in the 



field of view A is so small that the expected number of 
pairs would be unacceptably small. Combining data from 
M statistically independent snapshots reduces the asso- 
ciated error in g(r) to A^ s 'g(r) = g(r)/(2nn 2 AMrdr). 

Errors due to uncertainties in particle location can 
be calculated as A^g(r) = 2d r g(r) Ax, where Ax is 
the error in locating a single particle's centroid in each 
dimension. The radial derivative of g(r) can be com- 
puted numerically from the experimental data, which is 
binned to resolution dr Typically, Ax <C dr, so that 
A^g(r) < A^g(r). 

Even though the particles' out-of-plane excursions are 
small, they also contribute to errors in g(r) through pro- 
jection errors, especially near contact. Out-of-plane fluc- 
tuations Sh make particles appear to be closer than they 
actually are. The error in apparent particle separation 
falls off with separation as Ar — (6h) 2 /r. In practice, 
we combine this contribution in quadrature with the es- 
timated error due to inaccuracies in particle tracking, 
2 Ax, in computing A^g{r). 

Combining A^ m ^g{r) and A^g{r) in quadrature estab- 
lishes the range of possible values of g[r) for a given sam- 
ple, restricted only by the requirement that g(r) > 0. We 
compute trial pair potentials in both the HNC and PY 
approximations using both the upper and lower bounds 
on g{r) as inputs. The resulting lower and upper esti- 
mates on u(r) then are added in quadrature with the 
systematic error due to differences in HNC and PY re- 
sults to obtain estimates for the upper and lower error 
bars on u(r). Typical results appear in figure |3 and es- 
tablish that the minima reported in these data are indeed 
clearly resolved by our methods, even if the error bounds 
near contact are substantial. 



THERMODYNAMIC SELF-CONSISTENCY: 
CONFIGURATION AL TEMPERATURE 

Despite the care taken to estimate and eliminate 
sources of error in these measurements, using Eqs. © 
and J7J) to interpret experimental data might be criticized 
for its uncontrolled approximations: Eqs. © and Q can 
converge numerically to an answer even when applied well 
beyond their domain of validity. Assessing the bounds of 
this domain can be problematic if the form of the pair 
potential is not known a priori. Applying liquid struc- 
ture theory to experimental data also requires the as- 
sumption of pairwise additivity. Nonadditivity, however, 
would have no obvious signature in the results. Other 
unintended processes such as noncquilibrium hydrody- 
namic coupling also can yield reasonable-looking results 
that could be mistaken for an equilibrium pair interac- 
tion [sOl • Consequently, the appearance of qualitatively 
new features in any particular measurement of u(r) could 
signal a failure in the method. For this reason, most pub- 
lished accounts have relied upon comparisons among sev- 
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eral related systems to bolster their conclusions regard- 
ing trends in confinement-mediated interactions. These 
comparisons are themselves subject to question because 
the ultraclean chemical environments required for these 
studies are difficult to alter in a predictable manner. 

To address all such concerns, we have introduced |3| 
methods to assess whether or not a trial pair potential 
describes a system's interactions in a thermodynamically 
self-consistent manner. Our approach is based on the re- 
cently introduced notion of a configurational temperature, 
which has found widespread applications in simulations 
[3ll l32| , but has not previously been applied to experi- 
mental data 1241. 
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FIG. 3: Various consequences of the generalized tempera- 
ture definition for selected choices of the arbitrary vector field 
B(T). 

The temperature of an equilibrium ensemble of parti- 
cles is defined conventionally in terms of the particles' 
mean kinetic energy, without regard for their instanta- 
neous positions. In 1997, Rugh pointed out that the tem- 
perature also can be expressed as ensemble averages over 
geometrical and dynamical quantities [33j . This notion 
is expressed more generally [3, as 



(Vft(r)-B(r)) 
(v-B(r)) ' 



(9) 



where angle brackets indicate an ensemble average, T = 
{qi, . . . , <?3at,Pi, . . • ,P3n} is the instantaneous set of 3iV 
generalized coordinates qj and their conjugate momenta 
Pj for an iV-particle system, H(T) = Y^=i Pj I (2m) + 
V({qj}) is the Hamiltonian associated with the conser- 
vative iV-particle potential V({qj}), and B(T) is an ar- 
bitrary vector field selected so that both the numera- 
tor and denominator of equation are finite and the 
numerator grows more slowly than e N in the thermo- 
dynamic limit. Choosing B(T) = {0, . . . , 0,pi, ■ ■ -P3n} 
yields the familiar equipartition theorem. Choosing in- 
stead B(r) = — VV r ({<7i}) yields a formally equivalent 
result, 



^config 



(iwi 2 ) 



(10) 



which depends only on the particles' instantaneous con- 
figuration, and not on their momenta. 
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FIG. 4: (a) Estimating the range over which interactions 
affect the configurational temperature for the silica data at 
H = 9 fim. (b) Finite-size scaling of several variants of the 
configurational temperature for the silica data at H = 9 fim. 
Solid curves are fits to third-order polynomials in 1/N show- 
ing extrapolations to the thermodynamic limit. 

Directly applying equation (|10fl requires the full N- 
particle free energy, which is rarely available. Simplified 
forms emerge for systems satisfying certain conditions. 
For example, if V({qi}) is the linear superposition of pair 
potentials, u(r), then equation (flTTf) reduces to [sil ]. 



T N F'f 



(11) 



where Fj = — ^i u ( r ij) i s th 6 total force on particle 
i due to its interactions with other particles, V» is the 
gradient with respect to the i-th particle's position, r^, 
and ry = |rj — r,-| is the center-to-center separation be- 
tween particles i and j. The temperature is reflected in 
the instantaneous distribution of forces because objects 
explore more of their potential energy landscape as the 
temperature increases. 

Equation (|llf) may be generalized into a hierarchy of 
hyperconfigurational temperatures by choosing B(T) = 

{/;!: 



T N F 



8+1 



V,; • F, 



(12) 



for s > 0. These higher moments are more sensitive 
to the input potential's detailed structure than r con p = 
TP' . They also can be applied to three-dimensional sys- 
tems with long-ranged 1/r potentials, for which T con F is 
ill-defined. Equations © through Ijl2|l apply only in the 
thermodynamic limit, with errors of 0(1/N). 

For systems with short-ranged potentials, dropping ad- 



6 



ditional terms of 0(1/N) from equation © yields |34j |: 

/ y N f 2 \ 
fcsT COI1 i = -( vrw' = v,,V, / > anrf ( 13 ) 

/-■^....•j - ( >■'',; ■) , (14) 



\ £f =1 if / ' 

the second of which was proposed in reference [24| . These 
definitions' different dependences on sample size N are 

( s) 

useful for comparison with ' . 

We apply the configurational temperature formalism 
to our colloidal monolayers by using the measured par- 
ticle locations p(r,t) and extracted pair potential u(r) 
as inputs to the various definitions. Provided that the 
conditions for the configurational temperatures' deriva- 
tion are met, then all variants will yield results consistent 
with each other and with the (known) temperature T of 
the heat bath. In particular, consistent results emerge 
only if the system is in local thermodynamic equilibrium, 
if its interactions are indeed pairwise additive, and if the 
measured pair potential u(r) accurately reflects those in- 
teractions. 

In practice, each snapshot of a monolayer's configu- 
ration constitutes a measurement of its configurational 
temperature. Particles near the edge of the field of view, 
however, may have strongly interacting neighbors just 
out of the field of view whose contributions to their net 
force would be overlooked. Including these apparently 
unbalanced forces would grossly distort estimates of the 
configurational temperature. To avoid this, we calcu- 
late force distributions only for particles whose relevant 
neighbors all lie within the field of view. Such particles 
lie no closer than the interaction's range R to the edge 
of the field of view. We estimate R from u(r) and g(r) 
by computing 



T(r) 
T 



„ „ r , . |Vu(r)| 2 
2n(3-g(r) 1 V / , 
a V u[r) 



(15) 



an example of which is plotted in figure Ufa). This func- 
tion may be interpreted as the contribution to the con- 
figurational temperature due to particles separated by 
distance r. Quite clearly, pairs with r > R contribute 
little if at all to the configurational temperature. 

This necessary step further reduces the number N of 
particles in the field of view. This is problematic because 
all of the temperature definitions involve approximations 
of 0(1/N). Adopting a standard technique from sim- 
ulation studies, we deliberately subsample the available 
data, recalculate the configurational temperature on the 
restricted data set, and extrapolate to the large N limit 
by fitting the result to a polynomial in 1/N. Typical 
results appear in figure Hfb). Even though the different 
definitions have substantially different dependences on 
sample size, they all extrapolate to the thermodynamic 
temperature in the thermodynamic limit. 

This result turns out to be reassuringly sensitive to 
details of the pair potential. The small residual scatter 



in the experimental u(r) is greatly magnified in calculat- 
ing the configurational temperature, particularly for the 
higher-order hyperconfigurational temperatures. Conse- 
quently, the data in figure [5] were fit to a fifth-order 
polynomial whose coefficients were used in calculating 
figure Varying the pair potential by as little as one 
percent in the region of the core repulsion increases the 
apparent configurational temperature by more than ten 
percent. Simply truncating the attractive minimum in 
u(r) to mimic a purely repulsive potential leads to a fifty 
percent increase, or an error of 150°C. 

The successful collapse of the configurational and hy- 
perconfigurational temperatures to the thermodynamic 
temperature constitutes a set of stringent internal self- 
consistency tests for the accuracy of the measured pair 
potential and its correct interpretation. When combined 
with the considerations from the previous sections, we 
can improve the estimated resolution of our pair poten- 
tial to roughly 1/20 k B T. 

The observed confinement-induced attractions there- 
fore should be considered a real, pairwise additive con- 
tribution to the monolayers' free energy, at least in this 
range of ionic strength and areal density. Attractions of 
0.3 k B T are not strong enough to induce phase separation 
at such low areal densities, moreover. This is consistent 
with the assumption underlying the liquid structure for- 
malism that the system is in a single homogeneous phase. 



THE ROLE OF CONFINEMENT 

We next investigate the role of geometric confinement 
in inducing like-charge attractions in sedimented silica 
monolayers. Figure shows data from five different 
monolayers of silica spheres (a — 1.58 /im) in slit pores 
ranging in depth from H — 195 /im down to H — 3.2 /jm. 
FigureEIb) shows the associated configurational temper- 
atures. For all inter- wall separations, the monolayer is 
sedimented at roughly h = 900 nm, with the only obvi- 
ous difference being the inter-wall spacing. Well-resolved 
attractive minima are evident for plate separations as 
large as H — 30 /im. 

This observation contrasts with measurements on more 
highly charged polystyrene sulfate spheres, for which 
anomalous attractions appear only when the spheres are 
rigidly confined to the midplane, at separations no larger 
than H = 4a || . This difference may be due to the silica 
spheres' proximity to the lower wall. Why then would 
attractions not be evident at H = 200 urn [lfiL El E3|? 
More to the point, why would a second wall at H = 20<r 
make a difference? Trends in figure Ufa) suggest an ex- 
planation. 

One prominent feature of these data sets is that the ap- 
parent range of the core repulsion moves monotonically to 
smaller r as the inter-wall separation decreases. This dif- 
fers with the results of optical tweezer measurements on 
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FIG. 5: (a) Measured pair potentials for monolayers of sil- 
ica spheres 1.58 fj,m in diameter sedimented to a height 
h = 900 ± 100 nm above a glass surface, for a variety of inter- 
wall separations, H . The solid curves are fits to equation l|16|l 
with parameters tabulated in Table[I] The inset schematically 
represents the space charge model for confinement-induced 
attractions, (b) Configurational temperatures for each of the 
interaction measurements. 



polystyrene spheres, in which the depth of the attractive 
minimum varies with H, but not the range of the repul- 
sion . It is tempting to ascribe the trend in our silica 
data to a decrease in the effective Debye-Hiickel screen- 
ing length as the ratio of surface area to volume increases 
and diffusive contact with the ion exchange reservoirs di- 
minishes. If this were the case, however, we would expect 
the slope of u(r) near contact to decrease monotonically 
also. Instead, there is no discernible trend, presumably 
because the base ionic strength varies randomly from run 
to run. 



Referring to the DLVO result in equation JHJ for guid- 
ance, it would appear that the spheres' effective charge 
Z is the only other parameter that might be free to vary. 
Such variation is consistent, at least qualitatively, with 
predictions of charge renormalization theory 01 f° r silica 
spheres near charged silica surfaces. It would also explain 
the different behavior of polystyrene spheres whose more 
acidic surface groups are not so susceptible to charge reg- 
ulation by nearby surfaces |l7| . It leaves open the ques- 
tion, however, of why an attraction appears at all. 



SPECULATION: SPACE-CHARGE MEDIATED 
ATTRACTIONS 

A variety of mechanisms beyond Poisson-Boltzmann 
mean field theory have been proposed for confinement- 
induced attractions among like-charged colloid. These 
include attempts to compute London-like attractions due 
to fluctuations in the distribution of simple ions around 
the large spheres [3(| and density functional analysis of 
high-order correlations in the distribution of large and 
small ions 0,E3- The few that appear to reproduce ex- 
perimental observations have proved controver- 
sial [13, and none of the more widely accepted calcu- 
lations predicts an attraction of the range and strength 
observed experimentally, particularly if the simple ions 
are monovalent. Nor have computer simulations yet been 
able to address the regime of large charge asymmetry that 
appears to be necessary for this effect. Other approaches, 
however, may shed light on these anomalous interactions. 

The Kirkwood-Poirer formulation of electrolyte struc- 
ture for example, suggests that the correlations 
between macroions and simple ions can become non- 
monotonic in the strongly coupled regime. Hastings 
subsequently pointed out that these correlations in the 
simple ion distribution would lead to local violations of 
electroneutrality in regions between macroions [42j . and 
that the resulting effective interaction between macroions 
would include an attractive component. This result 
parallels the more recent thermodynamically consis- 
tent liquid structure calculation by Carbajal-Tinoco and 
Gonzalez- Mozuelos [ii^ . 

If we hypothesize that the distribution of counterions 
extending away from a charged surface also develops re- 
gions of space charge when modulated by nearby spheres, 
then the effective inter-sphere interaction should include 
a term accounting for sphere-space charge-sphere bridg- 
ing. In the absence of a theory for the actual simple ion 
distribution, we model the space charge's influence as 
the screened coulomb interaction between the spheres' 
effective charges and a point charge of valence q centered 
between them: 



/3u(r) = Z 2 \i 



( exp(Ka)^ exp(— kt) 
V 1 + 



exp{na / 



exp (- 



1 + na 
(16) 

Fitting the data in figure Ufa) to this form yields re- 
markably good agreement, with fitting parameters tab- 
ulated in Table IJ The screening lengths in all cases are 
consistent with the expected micromolar ionic strengths 
of our apparatus. The spheres' effective charge number 
appears to decrease systematically with wall separation 
in a manner at least qualitatively consistent with charge 
regulation theory 0] . Most tellingly, the effective space 
charge number is consistent with q = 10 at all separa- 
tions. If this model is to be taken seriously, this result 
suggests that the sedimented silica spheres are indeed 
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influenced by the nearby wall's counterion distribution, 
and that the resulting attraction is evident only when 
the core electrostatic repulsion is not too strong. Re- 
ducing the spheres' effective charge exposes the nascent 
attraction in this scenario. For the more highly charged 
polystyrene spheres, reducing the wall separation has lit- 
tle effect on the spheres' effective charge or the screening 
length, but increases the concentration of counterions be- 
tween the spheres. 



H (fj,m) 


Z 


k 1 (fim) 


q 


195 


7000 ± 400 


200 ± 20 


3±3 


30 


2500 ± 150 


160 ± 20 


10 ±2 


18 


2400 ± 150 


140 ± 20 


13 ±3 


9 


800 ± 100 


150 ± 25 


13 ±4 


3.2 


800 ± 100 


150 ± 20 


11 ±2 



TABLE I: Interaction parameters obtained from fits to the 
space charge model. 

This simple space-charge model appears to account for 
the available observations of like-charge attractions be- 
tween confined charge-stabilized spheres. Its interpre- 
tation points toward a correlation-based explanation for 
the effect, albeit of an extraordinary range. The mea- 
surements described in the present work should help to 
eliminate any remaining concerns regarding the validity 
and accuracy of the larger body of measurements in the 
literature, and their interpretation. The thermodynami- 
cally self-consistent measurement protocol we introduce 
should also find applications in the broader context of 
experimental soft matter research. 
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